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Abstract 

We derive covariant baryon wave functions for arbitrary Lorentz boosts. Model- 
ing baryons as quark-diquark systems, we reduce their manifestly covariant Bethe- 
Salpeter equation to a covariant 3-dimensional form by projecting on the relative 
quark-diquark energy. Guided by a phenomenological multigluon exchange repre- 
sentation of a covariant confining kernel, we derive for practical applications explicit 
solutions for harmonic confinement and for the MIT Bag Model. We briefiy com- 
ment on the interplay of boosts and center-of-mass corrections in relativistic quark 
models. 
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One of the major goals of modern electron and hadron accelerators is the investigation 
of the internal structure of hadrons, in particular of baryons: detailed information is 
extracted from scattering experiments at large momentum transfers of typically 1 GeV/c 
and beyond. The corresponding form factors map out the various internal (generalized) 
charge distributions and provide stringent information on the underlying quark and gluon 
degrees of freedom. Presently various experiments are ongoing with electron or photon 
beams at MAMI,ELSA, MIT, JLAB und DESY (1) and with proton beams at COSY and 
CELSIUS (2) and other labs. 

Practical calculations of form factors suffer in general from the pertinent problem of 
center-of-mass (CM) corrections for the many-body problem and from drastic effects from 
Lorentz contraction at increasing momentum transfers. While for the CM corrections var- 
ious recipes have been developed and applied in practical calculations (3-6), less progress 
has been achieved in the formulation of covariant baryon wave functions suitable for prac- 
tical calculations (7-24). A possible alternative, the evaluation of formfactors on the light 
cone, where Lorentz boosts are completely kinematical, has so far entered only selectively 
in practical applications at low scattering energies, beyond that such an approach suffers 
from other decreases, such as the loss of strict rotational invariance (25). As in general 
the construction of boosted, Lorentz contracted wavefunctions is nearly as complicated 
as the solution of the full problem, in most practical applications ad hoc and purely 
kinematical prescriptions for the rescaling of the coordinate along the direction of the 
momentum transfer are applied (examples are given ref. (26-27)). Thus, specific ques- 
tions, as the dependence of Lorentz corrections on the confining kernel in quark models, 
are not addressed. In addition, to hopefully minimize the influence of Lorentz contrac- 
tions formfactors are in general evaluated in the Breit frame, though experimentally they 
are measured in the lab system. 

In this note we formulate an economical model for covariant baryon wave functions, which 
leads to results suitable for practical applications. As it our main goal to end up with 
analytical formulae, we model the baryon - in the following we use the word proton, 
though our approach is fairly general - as a quark-diquark system and restict ourselves, 
without any loss of generality, to spin-isospin scalar diquarks (28). 

Our starting point is the manifestly covariant 4-dimensional Bethe-Salpeter equation (29) 

T = K GT and ^ = GT (1) 

with the vertex function and the Bethe-Salpeter amplitude T and respectively, and the 
interaction kernel K. In the two-body Greens function for the quark with mass m and the 
diquark with mass m* we fix the relative energy dependence from the covariant projection 
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on the diquark (30) 

which results up to (^g) in the single particle Dirac equation for the quark for systems 
with arbitrary overall 4-momenta Q = (E(P) = \J P 2 + M 2 , 0, 0, P) 

[eJF) ( e + ¿^"^)" (aq + /3m) ) ^ (g ' q) = 



^-|^(g, q ,k)^(g,k)ák (3) 



Without any details we add a brief comment on the CM corrections in our model: evi- 
dently there is a direct coupling between the internal and external momenta q and P, or 
equivalently, between boosts and the CM motion. In the rest system the leading center- 
of-mass corrections are absorbed for e = m + e&, where e& is the binding energy of the 
quark in 

+ ^ - K(r)) y?(r) = (4) 

with the reduced mass 1/ ¡i = l/m+ l/(m + rrf) for an arbitrary quark potential V n (r) 
(a detailed discussion of CM corrections are presented in a separate paper). 

The decisive step for a practical model is the formulation of a covariant interaction kernel 
in eq. (4). As the dynamics of the quark - quark interactions, particularly the microscopic 
nature of the confinement, lacks an understanding on the fundamental level of QCD, 
all current models in practical calculations rely on phenomenological formulations of the 
interaction kernel. Being unable to do better, we proceed here along similar lines: we 
assume that the interaction kernel can be presented as a superposition of appropriately 
weigthed gluon exchange contributions; quantitative parameters can be extracted in com- 
parison with studies to baryon spectroscopy, decay rates or form factors (31). Thus we 
start from the general kernel 

for arbitrary powers of n (which reflect different parametrizations of the confining kernel) 
(eq. (5) contains the linear confinement in the Cornell potential (32)). Upon projecting 
out the relative energy dependence this yields the covariant, 3-dimensional kernel 

(d \ n 1 
~A~2 \2ÍD\ 2 i 2~ 2 ~ ( 6 ) 

with the "quenching parameter" 

X(P) = M/VM 2 + P 2 = M/E(P) (7) 
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where we introduced the mass scale ¡i (to regularize the Fourier transform to coordinate 
space). Already simple power counting signals, that a kernel with the power n leads 
to confinement with ~ r 2n . Upon performing the corresponding Fourier transform to 
coordinate space and performing the limit /i^Owe find 

K n (P, q) -> (1 + (3) /2 V n {\/ (z/ A(P)) 2 + p 2 ) (8) 

where we introduced for convenience the particular Dirac structure of the kernel to facil- 
itate the evaluation of the resulting Dirac equation. Eliminating the small component in 
eq. (3) with the kernel from eq. (6) and upon dropping CM corrections and e 2 terms for 
compactness, we end up with the Schroedinger type equation for the large component of 
the Dirac equation 

p 

(2me b - \ 2 (q z - jj™) 2 ) - <¿ - 

- V n ^(z/\(P)) 2 + p 2 )/i?) u(z, p) = (9) 

(with the typical length scale R; in the rest system the equation above reduces to the 
standard spherical Schrdinger type equation for a particle with mass m). The final equa- 
tion defines with its connection to the small component by a simple differentation the full 
relativistic covariant quark - diquark wave function for arbitrary Lorentz systems. In the 
equation above we see the shortcoming from the phenomenological nature of the inter- 
action kernel: we absorb the explicit e and P dependence of the kernel in the definition 
of the energy scale V n for the confining force; including an explicit P dependence in V n 
would require a detailed knowledge of its microscopic origin. 

Approximate or numerical solutions for eq. (9) can be obtained for different confining 
szenarios (a more detailed investigation, such as also of the popular linear (heavy quark) 
confinement (33), is presented elsewhere). Here we enter only briefly into two szenarios, 
which allow a rigorous analytic solution for arbitrary systems: i. e. harmonic confinement 
and bag models in the limit n — » oo in eq(6). 

- Harmonic confinement: 
With the harmonic kernel defined as (34) 

K(P,q,k) = -12/tt iim ((á/d/i 2 ) 2 (/i/2 + (d/dp 2 )^ /3) 

- - ((l/\(P)) 2 (d/dq z ) 2 + (d/dq 2 ± )) 5(q z - k z )5(q ± - k ± ) , (10) 

the solutions for arbitrary excitations of the baryon are easily obtained in momentum 
space. After a redefinition of the longitudinal momentum and upon separating the 
longitudinal and the perpendicular component, the general solution is given by a 
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product of confTuent hypergeometric functions (35). Here we focus only on the 
nucleon as the quark - diquark ground state and obtain explicitly 



u(q z , q± )=Ne-^ (A 2 fe-t^)W 



(11) 



with the oscillator parameter a 2 = -^=, with V c oc 1/R A being the confinement 
strength and with the ground state energy 

As expected the standard solution for the spherical harmonic oscillator is recovered 
in the rest system, i.e. for P=0 and A(0) = l.As the characteristic result we find a 
quenching of the effective P-dependent with size parameter 

a\Pf = (X(P)af = a 2 (13) 

which leads to Lorentz quenching in coordinate space along the z - axis and thus to a 
significant increase of the longitudinal high momentum components with increasing 
P (Fig.l(a,b)); 

Bag Model: 

As mentioned above we generate the Bag from the transition n — > oo in the power 
of gluon-exchange kernel. As we are unable to present an analytical solution for 
arbitrary n (a closed solution for the z-component exists only in the limit of vanishing 
binding e& — > (35)) we first perform the limit n — > oo and then solve the equation 

P \ 
e b + — q z + m- (aq + (3m)J u(z,p)=0 (14) 

with the standard MIT boundary condition for the large and small components at 
z = X(P)R for the bag radius R. For the large component the ground state solution 
can be represented as 

u(z, p) = N cos(y z) J (k±p) (15) 

where the A dependence of the z-component again reflects the quenching of the bag 
(The extension to excited baryon states again is straightforward). The quenching 
fo the bag along the boost momentum is also reflected in the boundary condition 
z = XR for p = 0, which for the deformed bag can be solved only numerically (36). 
Characteristic results for 3 different boost momenta are presented for the large and 
small component of the bag ground state solution in Fig. 2. 
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Comparing our findings with current more phenomenological recipes we find that a 
general and simple extension of the parametrization of the spherical wave functions 
and momentum distributions in the rest system to a boosted system, by rescaling 
the size parameter of the system, but keeping otherwise the spherical character of 
the solutions, is certainly very unsatisfactory and breaks down completely for boost 
momenta of typically P/M > 1. Only for very small boost momenta P simple 
approximations, such as 

u(z, p,a) = u (r,a/ (VZ A(P))) and (16a) 

u(z,p,R) = exp(-(r/(-s/3 \(P)R)) 2 u(r, R) (16b) 

simulate very qualitatively Lorentz quenching of slowly moving systems. 

With increasing boost momenta the breaking of the spherical symmetry for the quenched 
bay leads for the ground (and all excited) state to the admixture of additional angular 
momenta, which drastically enhance the momentum spectrum of the ground state with 
increasing q. A characteristic ressult is shown in Fig. 3 for the d-wave admixture for 
different boost momenta. 

Summarizing our main findings in this note, we have formulated covariant wave functions 
and their transformation properties in an analytical quark - diquark model for the baryon 
and we find characteristic modifications from the baryon rest system to moving Lorentz- 
systems for different confining kernels. 

Our findings suggest possible extensions and basic shortcomings of the model. We feel 
that an extension of the model to mesons as qq systems, towards a more realistic quark- 
diquark description of baryons or to genuine 3-quark systems (together with a systematic 
inclusion of CM corrections) imposes only technical problems and is certainly feasible. 
Here we only mention that the quenching factor from eq.(7) is recovered in leading order 
for all current projections of the BS equation: as an example the Blankenbecler-Sugar 
reduction (30) for quarks with equal masses yields immediately 

G BBS (P,q)~5(q -P/Mq z ) (17) 

A more serious problem for confining kernels with a finite power in the interquark distance 
r is the precise formulation of Lorentz quenching for the kernels itself (in Bag models the 
dependence is absorbed in the boundary condition). Here the unsurmountable problem 
is our current lack in understanding the confining mechanism: it is not clear, how the full 
P dependence enters into the kernel (for an example compare ref.(37); however, different 
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szenarios may lead to quantitatively very different results for large P; for an example com- 
pare ref. (37)). We feel that a more realistic extension of present phenomenological quark 
models undoubtedly requires a much deeper analytical understanding of confinement. 
Here significant progress in various directions has been achieved recently, to mention only 
the modelling of confinement of QCD in the Coulomb gauge (38) or the extension of 
conecpt of instantons to merons as solutions of the classical QCD equations (39). 
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Figure Captions 

Fig. 1: Z-dependence of the large component of the harmonic oscillator ground state in 
coordinate (a) and momentum space (b) for different boost momenta P = (thin 
line), M (middle line) and 2M (thick line). 

Fig. 2: (a) Quenching and boundary conditions for the bag along the boost momenta P=0, 
M, 2 M. The functions f(z) and g(z) denote the large and small components of the 
ground state wave function (for a bag radius R = 1 fm). 

Fig. 3: D-state admixture to the harmonic oscillator ground state. Compared are the s- 
wave distribution for P=0 with the d-state component for P=M and 2 M (middle 
and thick line, respectively). 
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